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ABSTRACT 

We propose a method to determine the cosmic mass density 12 from redshift-space distortions 
induced by large-scale flows in the presence of nonlinear clustering. Nonlinear structures in 
redshift space such as fingers of God can contaminate distortions from linear flows on scales 
as large as several times the small-scale pairwise velocity dispersion a v . Following Peacock & 
Dodds (1994), we work in the Fourier domain and propose a model to describe the anisotropy in 
the redshift-space power spectrum; tests with high-resolution numerical data demonstrate that 
the model is robust for both mass and biased galaxy halos on translinear scales and above. On 
the basis of this model, we propose an estimator of the linear growth parameter (3 = 
where b measures bias, derived from sampling functions which are tuned to eliminate distortions 
from nonlinear clustering. The measure is tested on the numerical data and found to recover 
the true value of /? to within ~ 10%. An analysis of IRAS 1.2Jy galaxies yields f3 = 0.8+0.3 
at a scale of 1,000 km/s which is close to optimal given the shot noise and the finite survey 
volume. This measurement is consistent with dynamical estimates of /? derived from both real- 
space and redshift-space information. The importance of the method presented here is that 
nonlinear clustering effects are removed to enable linear correlation anisotropy measurements on 
scales approaching the translinear regime. We discuss implications for analyses of forthcoming 
optical redshift surveys in which the dispersion is more than a factor of two greater than in the 
IRAS data. 

Subject headings: cosmology: theory — large-scale structure of universe — galaxies: 
clustering 


1 Introduction 

A redshift-space map of galaxies is distorted relative to the real-space galaxy distribution as a result of 
peculiar motions along an observer’s line of sight. These distortions generate an anisotropy in pairwise 
correlations which would not appear in the real-space galaxy distribution of a statistically homogeneous and 
isotropic universe. Kaiser (1987) quantified the correlation anisotropy that results from large-scale peculiar 
flows in terms of the power spectrum of galaxies using the linear theory of gravitational instability. He 
demonstrated that in the linear regime the power measured by a distant observer only depends on the angle 
between the wavevector and the observer’s line of sight, and on the dimensionless factor f3 = H°- 6 /6. Here, 
H is the cosmic mass density parameter and b is a function that differs from unity if the galaxies are a 
biased sample of the total mass. Hamilton (1992) transformed Kaiser’s result out of the Fourier domain to 
determine the redshift-space correlation function £ s in linear theory and proposed an estimator of (3 based 
on a spherical harmonic decomposition of £ s . Subsequently, a number of (3 measurements from linear flow 
distortions have been reported (e.g., Hamilton 1993, 1995; Bromley 1994; Fisher, Scharf & Lahav 1993; 
Fisher et al. 1994; and Nusser & Davis 1994). 

There are two major challenges for a statistical measurement of (3 from linear flows in a real galaxy 
redshift catalog. The first comes from the finite size of the catalog. The catalog must contain a fair 
sample of structure on a particular scale or else the correlation information will contain finite-sampling 
noise. A reasonable criterion is that sampling scales should be below ~ 10% of the characteristic size of 
survey volume. The second challenge is to ensure that the sampling scales are large enough so that signal 
from nonlinear clustering does not contaminate the linear fluctuation modes. Fingers of God can extend 
up to several thousand kilometers per second in optically selected galaxy surveys and clearly affect the 
redshift-space power on these scales. (Here redshift-space distances are usually given in terms km/s; where 
real-space lengths are needed, we use the conventional units of h -1 Mpc where h is the Hubble parameter in 
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units of 100 km/s.) Thus the signal in redshift space from linearly growing fluctuation modes on scales from 
~1,000 km/s (characteristic of the translinear regime below which linear theory breaks down) to ~4,000 km/s 
can be contaminated by strongly nonlinear features. Present-day catalogs are typically within a few times 
10,000 km/s in radial extent, leaving at best only a small range of scales on which purely linear modes can 
be measured with a high degree of statistical integrity. 

In this paper we address this problem of determining [3 from linear flows in the presence of nonlinear 
clustering. Our approach is motivated by a remarkable result from Peacock & Dodds (1994) concerning the 
redshift-space power spectrum P s . While £ s itself can be modeled given the distribution of pairwise velocities, 
there is evidence (Fisher et al. 1994; Warren 1995) that this distribution has complicated behavior on a 
critical range scales from ~1 10 ft -1 Mpc. In contrast, Peacock & Dodds (1994) found that nonlinear signal 
can be modeled simply and accurately in the Fourier domain on translinear scales and above. Specifically, 
the effects nonlinear clustering and linear flows in redshift space can be described by separable linear filters 
acting on the real-space power spectrum. The linear filtering hypothesis has been supported by tests with 
7V-body simulations, although only mass particles and not collapsed galaxy halos were considered (Gramann, 
Cen & Bahcall 1993; Tadros & Efstathiou 1996). Here, we first suggest a form of the filters, which amounts 
to constructing a model for the anisotropy in P s , and determine its validity in high-resolution numerical 
simulations using both mass particles and collapsed halos which are identified as galaxies. On the basis of 
this model we then suggest a method to extract (3 on scales of 1,000 km/s and above, even when the velocity 
dispersion is of comparable size. Finally, we apply our method to the IRAS 1.2Jy survey. The result is a 
new measure of (3 on translinear scales. 


2 Power in Redshift Space 

Following Peacock & Dodds (1994) we assume that the effects of line-of-sight peculiar motion on a galaxy 
distribution can be modeled by a linear filter function F that relates the real-space and redshift-space power 
spectra: 

P s (k,ii) = F(k, n)P R (k) , (1) 

where /i is the angle cosine of the observer’s line of sight and a wavevector of amplitude k. For definiteness 
we assume that the power is measured in a distant compact survey so that /i is constant for all points in the 
survey. 

The form of F can be motivated from its limiting behavior at small and large scales. At small scales, the 
peculiar motions are dominated by isotropic dispersion, and the redshift-space correlation function obtains 
from a simple convolution of the real-space function and the distribution function f v of pairwise radial 
velocities. According to evidence from simulations (Zurek et al. 1994) and observations (Marzke et al. 1995), 
f v is accurately described by an isotropic exponential function at small pair separation, 

fv ~ exp(-v / 2|u| /o v ) , (2) 

where cr v is the pairwise radial velocity dispersion. In the Fourier domain, the redshift-space power spectrum 
is thus related to the real-space power spectrum by a linear filter, F D , the Fourier transform of f v taken as 
a function of line-of-sight pairwise velocity: 


F d [k, n) 
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( 3 ) 


Turning to the large-scale regime, Kaiser’s (1987) formula for a filter F L which accounts for the anisotropy 
in redshift-space power from linear flows is 


F L =(l+f3^ 2 ) 2 . 


( 4 ) 


Peacock & Dodds (1994) found that an effective model for the desired filter F is obtained from a product 
of F d and F L . Following their ansatz we write equation ®) as 
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In committing to this model (eq. ||) we limit our analysis to scales above the translinear regime. There is 
evidence (Brainerd et al. 1996) that Kaiser’s formula (eq. ||) for the anisotropy in the redshift-space power 
spectrum breaks down at translinear scales. This breakdown can be quantified by including a “softening 
parameter” of roughly translinear scale into equation (||) which drives the filter to unity in the strongly 
nonlinear regime (see equation [5] in Brainerd et al. 1996). Theoretical insight into this effect is given by 
Fisher & Nusser (1996). 

The model in equation (]s|) also may be limited by the assumption that the small-scale velocities are 
random and isotropic and exhibit none of the coherent flows that are generated by the linear or quasilinear 
dynamics. Numerical simulations support this assumption on megaparsec scales and below. A naive linear 
theory extrapolation of the pairwise velocity dispersion by way of the cosmological pair conservation equation 
(Peebles 1980 §71) also suggests that the linear flow contribution to o v at 1 h~ x Mpc is small for most 
plausible cosmological scenarios. In the cold dark matter model discussed below, linear flows would contribute 
(in quadrature) only ~ 100 km/s to the total of cr v ~ 1,100, which amounts to less than a percent in 
<r„. However the situation is more complicated on translinear scales where both random dispersion and 
significant coherent flows exist. Fortunately, the random motion contribution to the anisotropy of redshift- 
space fluctuation modes above the translinear regime comes primarily from strongly nonlinear scales where 
the velocity dispersion and spatial pairwise correlations are strongest and where the velocity distribution 
function takes a particularly simple form. 

We note that the above expression differs from the result of Peacock & Dodds in the choice of the model 
for the small-scale velocity distribution function. Their work was based on a Gaussian model, however, as 
they and others (Cole, Fisher & Weinberg 1995; Brainerd et al. 1996) have noted, the difference between a 
Gaussian distribution and the exponential distribution is of little consequence in the quantification of redshift- 
space distortions. Nonetheless we advocate the exponential model because of its successes in describing both 
real and simulated data. 

To determine the efficacy of the parameterization in equation (! we examine the redshift-space power in 
numerical simulations. We ran two 17 million particle high-resolution N-body simulations in periodic cubes 
of size 125 h~ x Mpc and 500 /i _1 Mpc, respectively, using a treecode with force smoothing of 0.01 h~ 1 Mpc. 
In both cases the primordial power spectrum was based on cold dark matter (CDM) normalized to erg = 0.74 
in an Q = 1, h = 0.5 universe. This cosmogony provides a reasonable match to the observed spatial clustering 
and pairwise velocity statistics of optically selected galaxies (Zurek et al. 1994; Marzke et al. 1995; Brainerd 
et al. 1996). We use the large-volume simulation to study the redshift distribution of mass on scales greater 
than 50 /i _1 Mpc; the small-volume run has sufficient mass resolution to identify halos as realistic galaxy 
candidates, hence it is used to study both the mass and galaxy distributions on smaller scales. 

The results of the N-body analysis are presented in Figure 1, showing the anisotropy in redshift-space 
power spectrum of the CDM mass (Fig. la) and the full catalog of halos (Fig. lb). We plot the power at 
constant k values as a function of the angle between the wavevector and the observer. Using the megaparsec 
separation value of cr Vl with values of 1,100 km/s and 850 km/s for mass and halos respectively, we find 
that the parameterization in equation (||) is successful on scales near and above the translinear regime. The 
breakdown in the model occurs on scales which are below o v , corresponding roughly to the characteristic 
size of the Fingers of God. 

Remarkable is the extent to which the dispersion signal affects the angular dependence of the redshift- 
space power. In the case of the CDM mass distribution, there is negligible anisotropy at 2x/k = 3, 300 km/s, 
three times cr v , because the large-scale enhancement in power along the line-of-sight is almost completely 
canceled by small-scale dispersion. 

Figure lc depicts results from a third dataset, a cool subset of halos obtained by eliminating all objects 
which are within 3 Abell radii of the highest density peaks and which have peculiar velocities greater 
than 750 km/s. Although produced in a rather contrived manner, this culled dataset is suggestive of the 
IRAS 1.2Jy galaxies in that they both have low velocity dispersions compared with optical galaxies and have 
relatively small numbers of objects in the cores of rich optical clusters. Indeed, the resulting catalog has o v 
of 384 km/s and spatial clustering properties below 100 h~ 1 Mpc which are consistent with the IRAS 1.2Jy 
data (Brainerd et al. 1996). Because of an antibias in the culled halo distribution relative to mass, (3 is about 
~ 10% above unity on translinear to linear scales with weak scale dependence. Given the great difference 
in small-scale clustering statistics between this dataset and the total mass distribution, it is significant that 
equation (! still gives a good parameterization of the anisotropy in redshift-space power on scales as small 
as 1000 km/s. 
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Fig. 1. -The two-dimensional redshift-space power spectrum. The power is shown as a function of the angle between 
the observer’s line of sight and the wavevector for several different wave magnitudes k. The plots show the power for 
(a) the mass distribution, (b) the galaxy halos, and (c) the culled halo subset with low small-scale velocity dispersion. 
The solid lines show the estimate of P from the simulation data, while the dashed curves are the predicted power 
using the multiplicative filter model in equation (5). A (3 value of unity was used to generate the model curves in 
all three plots, while the <r v values were 1,100 km/s, 850 km/s, and 400 km/s, for the mass, halos, and cool halos, 
respectively. 
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3 Dispersion-Corrected Sampling Functions 

With a working model for the anisotropy of redshift-space power in the presence of linear flows and small- 
scale dispersion, we now focus on measuring [3 from the linear contribution. The strategy is to determined 
the small-scale o v in a redshift-space density field, remove the effects of the dispersion signal as modeled by 
equation (^]), and then to estimate the anisotropy in Fourier modes on scales >10 /i _1 Mpc. In an infinite 
or periodic dataset one could simply remove the dispersion signal by a mathematically trivial deconvolution. 
For example, the redshift-space distribution might be transformed into the Fourier domain and the result 
multiplied by the inverse square-root of F D . However, in practice the finite size of a survey complicates 
the calculation of the Fourier transform especially at large wavelengths (e.g., Fisher et al. 1993). Here we 
circumvent the problem of the Fourier transform and deconvolution of the density field by making use of 
sampling functions which are tuned to be orthogonal to the dispersion signal. We consider only the variances 
in such samples and hence we extract linear flow information without directly deconvolving the redshift-space 
density field. 

A sampling function s gives a sample value S at some randomly chosen origin by projection onto a density 
field p 

S = Jd 3 x s(x)p(x)/p , (6) 

where p is the mean density. It is straightfoward to show that the variance in a set of many samples taken 
at random locations in a density field is related to the power spectrum P by 

ct 2 = Jd 3 kP(k)\s{k)\ 2 . (7) 

Of course in practical applications the power is inferred from a discrete galaxy distribution which contains 
shot noise. In this case, equation (iTj) must be corrected by subtracting the shot noise contribution, 

= [d 3 k\s(k)\ 2 , (8) 

n J 


where n is the mean density of galaxies. 

Anticipating the case when P in equation ( 0 ) is set to P s , we choose s to have a form such that in the 
Fourier domain 


(k) = s L (k,p) 


( 9 ) 


where s L couples only to the signal from linear flows. We assume that the samples are taken at a large 
distance from the observer so that p is constant. Clearly, the effect of the denominator in equation (|ii|) is to 
perform the deconvolution of dispersion signal at the level of random sampling. 

We emphasize that in the method discussed here we evaluate only sample variances; we never directly 
perform a deconvolution of the redshift-space density field to remove shot noise and the effects of small-scale 
dispersion. Yet the filter given in equation (jf|) might be used to do just that. Indeed, deconvolution with 
this filter and the inverse of the linear flow filter in equation (|J) would enable the real-space density field to 
be reconstructed from redshift-space data even on linear scales which are affected by small-scale dispersion. 
This would advance known reconstruction techniques (e.g., Yahil et al. 1994; Nusser & Davis 1995; Tegmark 
& Bromley 1996) by introducing some correction for nonlinear dispersion. However, the filtering must be 
done with great attention to noise and finite boundary effects. For example, optimal (Wiener) filtering is 
necessary to avoid the pitfalls of deconvolution in the presence of noise (cf. Press et al. 1992, section 13.3). 

With the nonlinear effects removed, the remaining problem is to determine a measure of (3 from the 
linear flow signal in the sample variances. Bromley (1994) proposed an estimator that follows from taking 
the derivative of some radially symmetric function R along the direction of a unit vector n at fixed orientation 
relative to the distant observer: 


s L (x) = n ■ X7R(x) and s L (k) = in ■ k R(k) . ( 10 ) 

Here R is specified to have the form of a ID Gaussian, 

R(x) = exp(— x 2 /2\ 2 ) , 


( 11 ) 
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Fig. 2.—The recovered (3 from the ratio-of-sample-variances method. In the plots the upper curves were obtained 
using a a v of 1,000 km/s, 850 km/s, and 425 km/s, respectively, and the lower curves have no correction for dispersion 
( <7 V = 0). The results for the mass distribution are shown in (a) while (b) contains the results for the full halo catalog 
(solid curves) and the culled halo subset (dashed curves). 



where A defines the scale of the sampling function. The effects of linear flows in redshift space can be 
isolated by considering two different sampling functions constructed from equations (p[) through P), with 
n directed parallel and perpendicular to the observer’s line of sight, respectively. The ratio of variances <7|| 
and <t_l from these two sampling functions is easily related to P using equation (7) and the redshift-space 
power in equation (j^): 


Q(f) = 4 


1 + f p + f/3 2 

i + 4 + i/? 2 ' 


( 12 ) 


The desired measure of p comes from inverting this expression. Note that the result does not require any 
knowledge the form of the real-space power spectrum as all dependence on P R cancels out in the ratio Q. 

We use the numerical data described in §2 to determine how well the ratio of sample variances in 
equation dH) is able to recover the correct value of p. Figure 2 illustrates the recovery from redshift 
space data over a range of sample function scales A using known values of the velocity dispersion a v on 
megaparsec scales from the three simulation catalogs. The method works well for the mass particles and 
the full halo catalog; in both cases the estimate of P was within 10% of its true value. In the culled halo 
subset the recovered P was consistently about 10% below the expected P value of ~ 1.1 when the a v value 
of ~ 384 km/s at 1 h Mpc pair separation was used in the sample functions. However, the number of halo 
pairs at a megaparsec is separation is relatively small, and with <j v itself rising fairly rapidly with separation 
at this scale, it is not surprising that a higher value of o v , 450 km/s, gives a better estimate of p. 

Figure 2 shows strong promise for determination of P with the ratio-of-variances method using tuned 
sample functions. The figure also illustrates the importance of taking into account small-scale velocity 
dispersion. Even on redshift-space scales approaching three times the pairwise velocity dispersion, we find 
that the recovered P is generally more than 50% below its true value. 


4 The IRAS 1.2 Jy Survey 

We now apply the sampling method to the IRAS 1.2Jy redshift catalog. This is a flux-limited sample of 
5,665 galaxies with sky coverage of 87.6% (Strauss et al. 1990; Yahil et al. 1991). The flux limit in the 
60/im waveband causes the density of galaxies in redshift space to fall off rapidly with distance; the number 
of galaxies in radial bins peaks at ~5,000 km/s, and the selection function <fi which describes the fall-off has 
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a value of less than 0.01 at 10,000 km/s. We take this latter distance to be the characteristic size of this 
survey. 

The sampling procedure used to generate Q from the IRAS galaxies is more complicated than in the case 
of the pristine numerical data. First, each galaxy must be weighted by l/</> to account for low-luminosity 
neighbors which were missed because of the flux limit. Here we use a selection function of the form given by 
Yahil et al. (1991, eq [11] therein) with a = 0.51, b = 1.84, and hr* = 5,440 km/s. 

The second complication in the analysis of the IRAS galaxies is that shot noise is significant on scales 
below ~1,000 km/s. If uncorrected it causes the sampling method to underestimate Q and (3. Here, we 
remove the shot noise contribution to the sample variance by equation (|^); its amplitude for each level of 
sampling scale is determined from artificial catalogs generated by a Monte Carlo method in accordance with 
the selection function of the IRAS galaxies. 

A third complication is incomplete sky coverage. So that we may sample with the tuned functions without 
the limitations of avoiding unobserved sky regions, we follow the strategy of Fisher, Scharf & Lahav (1994) 
who interpolated the IRAS galaxy density field across the zone of avoidance and other unsampled regions. 
The interpolation method is described in detail by Yahil et al. (1991) as part of a broader algorithm to 
reconstruct the real-space density field from the IRAS data in redshift space. Here, we include 705 artificial 
galaxies in real space kindly provided by M. Davis. Since the artificial galaxies represent a real-space density 
field, reconstructed under the assumption that (3 is unity, the line-of-sight mode amplitudes are guaranteed 
to be smaller than in the redshift-space IRAS catalog. This effect, along with the fact that any discontinuity 
with a boundary running along the observer’s line of sight will only add power to wavemodes in the plane 
of the sky, can drive the ratio Q to systematically lower values, causing (3 to be underestimated. 

There is one final difference between the analyses of the numerical data and of real galaxies. In the 
simulations the pairwise velocity dispersion a v at megaparsec scales, the essential ingredient of the dispersion 
filter (eq, j3j), is known precisely from the 6D phase space distribution of particles or halos. In contrast, a v 
must be gleaned from redshift-space data in the IRAS catalog, a procedure which introduces some uncertainty 
into a v . Published values (e.g., Fisher et al. 1993) and our own estimate based on the method of Davis 
& Peebles (1983) suggest that a v for the IRAS 1.2Jy sample is 320 km/s with 1-cr errors of ~ 40 km/s. If 
the IRAS galaxies are unbiased tracers of mass, then the numerical data suggest that this value should be 
appropriate for use in the dispersion filter (eq. |J). If the IRAS population is antibiased with respect to 
mass, as has been conjectured on the basis of comparison with optical surveys, then the use of this value may 
generate an underestimate of (3. This conclusion is made on the basis of similarity between our culled halo 
catalog and the IRAS galaxies which both have relatively low number densities in cluster cores (as compared 
to the full halo catalog or optically selected galaxies) and rapidly rising o v values on megaparsec scales. 

After implementing the above modifications to the sampling procedure, we analyzed the IRAS galaxies 
which lie inside a thick spherical shell of redshift space with an inner radius of 500 km/s and an outer radius 
chosen between 10,000-12,000 km/s. We estimate uncertainty by determining an approximate number of 
spatially independent samples that can be measured with a given sampling function and make repeated 
measurements of (3 with that number of samples. In Figure 3 we show the recovered (3 from the IRAS 
galaxies as a function of sampling scale. In general, the recovered value has a lower bound at the l-o level 
which is above 0.4. The mean of the recovered f3 on scales of 600-1,500 km/s is in the range of 0.6-0.9, 
depending on the outer radius of the sampled region. This sensitivity shows up mostly at the larger sampling 
scales where we must balance a desire to use as big a volume as possible against the increasing sparsity of 
galaxies at large distances. The most stable estimate comes from sampling scales around 1,000 km/s for 
which we find that (3 = 0.8+g'g at the 1 — a level. As noted above, the systematic error from zone-of-avoidance 
interpolation is likely to cause this value of (3 to be a slight underestimate of the true value. The similarity 
between the IRAS galaxies and our culled halo catalog in §3 mentioned above suggests that our choice of 
a v = 320 km/s may be too small, leading to a more significant underestimate of (3. Indeed, cr v in the IRAS 
data climbs to a value of 370 km/s at 2 h _1 Mpcwhich leads to a ^10% increase in (3 over the estimates 
shown in Figure 3. 


5 Conclusions 

Here we have modeled the anisotropy in the redshift-space power spectrum on the basis of Peacock & 
Dodds’ (1994) theoretical work and our own high-resolution numerical simulations. We demonstrated that 
the effects of nonlinear clustering on the power spectrum of both mass and galaxy halos can be described 
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Fig. 3.—The recovered (3 from the IRAS 1.2Jy survey. The upper curve comes from sampling with functions tuned 
to a v = 320 kms, and the lower curve corresponds to zero dispersion. The error bars show the 1 -<t uncertainties. 


by a simple multiplicative filter function. The filter-based model for the anisotropy is seen to hold on linear 
and translinear scales which are greater than the pairwise velocity dispersion a v at megaparsec. Having 
established the domain of validity, we then proposed a redshift-space statistic, the ratio of sample variances 
in equation & which is both impervious to the effects of nonlinear clustering and sensitive to linear flows 
and the linear growth parameter (3. 

This work should serve as a cautionary reminder of the extent to which the nonlinear distortions con¬ 
taminate signal from linear flows (e.g., Fig. 2). In a measurement of power, the magnitude of contamination 
can be determined from equation ©>; for example the nonlinear signal reduces the power along the line of 
sight by ~ 40% at scales 2ir/k equal to five times o v . In the case of the IRAS 1.2Jy survey where a v is 
~ 300 km/s, the presence of nonlinear distortion is important even at scales of 1,500 km/s. For optical 
surveys where a v ~ 800 km/s (e.g., Marzke et al. 1995), the extent of contamination exceeds 4, 000 km/s. 

The work presented here also gives the encouraging message that nonlinear clustering does not prevent 
the accurate extraction of (3 from linear flows in redshift data. Figure 2 provides the evidence based upon 
cosmological simulations of the greatest dynamical range published to date. 

The analysis of the IRAS 1.2Jy survey is intended to show that the method introduced here gives rea¬ 
sonable results. Our estimate of (3 = O. 8 +Q 3 is consistent at the 2-cr level with virtually every published 
value based on an analysis of the IRAS galaxies (e.g., Dekel et al. 1993; Fisher, Scharf & Lahav 1994; Cole, 
Fisher & Weinberg 1995; Hamilton 1995), although it is higher than most redshift-space measurements based 
exclusively on linear theory. Our procedure is most similar in spirit to that of Hamilton (1995) who obtained 
/3 = 0.69toug using a merged IRAS 1.2Jy and the optically selected QDOT catalog. Hamilton, working 
with harmonics of a smoothed power spectrum, modeled the nonlinear noise in the same way as we did, 
with an isotropic exponential distribution for pairwise velocities. We hope that we have provided convincing 
evidence in §2 that such a model is indeed robust. 

We note that our estimate of /? from the IRAS 1.2Jy sample is consistent with dynamical analyses that 
make use of real-space information as well as redshift data. For example, Dekel et al. (1993) find j3 = 1.28±0.3 
on scales of 1,500-4,000 km/s. While knowledge of both real-space and redshift distributions can yield a 
more accurate measurement of (3 , the real-space data are difficult to measure and therein lies the value in 
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extracting (3 directly from redshift space. Fisher et al. (1994) performed a pure redshift-space analysis of the 
IRAS 1.2Jy galaxies by modeling £ s and obtained P = 0.45+0^8 on scales of 1,000-1,500 km/s, well below 
the dynamical estimates at larger scales. However, the modeling of and its dependence on the uncertain 
behavior of the velocity distribution function is difficult. In contrast, our measure of P, based instead on a 
much simpler model in the Fourier domain, is consistent with the large-scale dynamical results. 

The ultimate hope is that our measure may be applied to forthcoming optical surveys which exhibit high 
o v values (e.g., Marzke et al. 1995). With the ability to accurately determine P down to translinear scales 
in spite of nonlinear clustering effects, our measure can provide high statistical significance in a correlation 
analysis of redshift-space data. 

We are grateful to Peter Quinn and the Mount Stromlo Observatory for hospitality at the Heron Island 
Workshop on Peculiar Velocities. This work was supported under the auspices of the U.S. Department of 
Energy and by the NASA HPCC program. The Cray Supercomputer used in this investigation was provided 
through funding from the NASA Offices of Space Sciences, Aeronautics, and Mission to Planet Earth. 
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